Unambiguous pure-state identification without classical knowledge 
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We study how to unambiguously identify a given quantum pure state with one of the two reference 
pure states when no classical knowledge on the reference states is given but a certain number of 
copies of each reference quantum state are presented. By unambiguous identification, we mean that 
we are not allowed to make a mistake but our measurement can produce an inconclusive result. 
Assuming the two reference states are independently distributed over the whole pure state space in 
a unitary invariant way, we determine the optimal mean success probability for an arbitrary number 
of copies of the reference states and a general dimension of the state space. It is explicitly shown that 
the obtained optimal mean success probability asymptotically approaches that of the unambiguous 
discrimination as the number of the copies of the reference states increases. 



I. INTRODUCTION 

In quantum mechanics one cannot perfectly clone an 
unknown state [lj, which makes the problem of distin- 
guishing quantum states nontrivial [3, Q • Imagine we 
are presented with an unknown quantum pure state p on 
a d-dimensional vector space C d . Let us assume that the 
input state p is guaranteed to be either one of two ref- 
erence states pi and p 2 , each being also a pure state on 
C d . Then how well can we identify the input state with 
one of the two reference states? 

We can consider two cases depending on what kind of 
information on the reference states is available. In the 
first case, it is assumed that we have complete classical 
knowledge on the two reference states p\ and pi. This 
is the standard setting of quantum-state discrimination, 
which was solved by Helstrom [2|. 

On the other hand, we can also consider the case where 
only a certain number (iV) of copies of p\ and p 2 are pre- 
sented, with no classical knowledge on them available 0. 
See also related works in the case of qubits In this 

case, we could obtain only limited classical information 
on the reference states, since the no-cloning theorem [lj 
does not allow us to increase the number of copies of 
the reference states. The best we can do is to perform a 
positive-operator-valued measure (POVM) measurement 
on the total state p ® pf N ® pf N and try to identify the 
input state p with one of the reference states p\ and p 2 ■ 
If the number of copies, TV, is infinite, the problem is 
reduced to quantum-state discrimination, since we could 
always obtain complete classical knowledge of a quantum 
state. In our previous paper p|, we called this problem 
"state identification" and determined the optimal mean 
identification probability for an arbitrary number (N) of 
copies of the reference states in a general dimension d. 

In the standard setting of the discrimination problem, 
we are allowed to make an error and are interested in the 
optimal strategy that makes the error probability mini- 
mum. On the other hand, an error is not allowed in the 
problem of unambiguous discrimination [3, la, |9( . Instead 
our measurement can produce one of three outputs 1,2, 
or 0. If the output is 1(2), we are certain that the in- 



put state p is pi(pz), and the output means that we 
do not know the identity of the input, which is called an 
inconclusive result. The optimal strategy is the one that 
minimizes the probability of the inconclusive result. 

We can also generalize the unambiguous discrimina- 
tion problem to the case in which a finite number (N) of 
copies of the reference states are presented without any 
classical information on them, which will be called the 
problem of unambiguous state identification in this pa- 
per. Bergou and Hillery studied this problem in the case 
of qubits (d = 2) when the number of copies of reference 
states, N = 1 10] . They called the optimal strategy a 
programmable state discriminator since the strategy is 
not "hard wired" but supplied by the reference states 
stored in registers in the machine. 

The problem of the quantum-state comparison [ill IT^ 
is related to the state identification problem. One's task 
here is to establish whether or not two quantum sys- 
tems have been prepared in the same state. The sym- 
metry under interchanging the systems is essential in the 
state comparison, since the combined system is symmet- 
ric when the two systems are in the same state, whereas 
it has no definite exchange symmetry otherwise. The ex- 
change symmetry plays a crucial role also in the state 
identification, but in a more involved way. 

In this paper we study the problem of unambiguous 
identification of pure states for an arbitrary number iV 
of copies of the reference states in a general dimension 
d. The two pure reference states p\ and p 2 are indepen- 
dently distributed over the whole state space on C d in a 
unitary invariant way. The input state p is assumed to be 
either one of the two reference states with the equal prob- 
abilities. We determine the optimal POVM and the opti- 
mal mean unambiguous identification probability, which 
can be explicitly shown to approach the mean unambigu- 
ous discrimination probability in the large- N limit. 



II. MEAN UNAMBIGUOUS PURE STATE 
DISCRIMINATION 

In this section we average the unambiguous discrim- 
ination probability, assuming that the two reference 
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states are independently distributed on C in a unitary- 
invariant way. This mean unambiguous discrimination 
probability will be later compared with the mean unam- 
biguous identification probability in the large- N limit. 

More precisely the unitary distribution of the reference 
states is specified in the following way. Expand a pure 
state as |0) — X)f=i c »l*) m terms of an orthonormal 
base {\i)} of C d . The distribution is then defined to 
be the one in which the 2d-component real vector {xi = 
Heci,yi = Imci} is uniformly distributed on the (2d— 1)- 
dimensional hypersphere of radius 1 with the integration 
measure given by 



dcdc+ = JlidxtdyJS foz? + yf) - l^j 



(1) 



Evidently the distribution does not depend on the choice 
of the orthonormal base {| i )}. 

The optimal success probability of unambiguous dis- 
crimination of two known pure states p\ = \ 0i ) ( 0i | and 
P2 = | 0i | on C d is given by Hill 



Pmax(Pl,P2) = 1 - 1(01 I 02 )|- 



(2) 



We calculate the mean unambiguous discrimination 
probability 



(3) 



where (• • •) means the average over pi and p2, which are 
independently distributed according to the unitary dis- 
tribution defined above. 

The average (\(4>i | 02 ) |) can be calculated in terms of 
an integration over a (2d — l)-dimensional hypersphere 
as 



(1(01 I02)|) = 



jdcdc + \c\ | 
Jdcdc+1 



(4) 



with the integration measure given by Eq. l[Tjl. 
We obtain 



Pmax(<2) = 1 



2 d ~\d- 1)! 
(2d-l)!! ' 



(5) 



which is certainly less than the mean discrimination prob- 
ability given in £|: 

& 1 + ^H<*ilfc>l 2 )H + ;^V (6) 



III. UNAMBIGUOUS PURE STATE 
IDENTIFICATION 

Suppose we are given an unknown pure state p on C d . 
We know that p is either one of the two reference states 
pi and pi on C d , with equal prior probabilities. Let 
us assume that we have no classical knowledge on the 



reference states, but a certain number (N) of copies of 
each state are available. What is the optimal strategy to 
unambiguously identify the input state with one of the 
reference states when the two reference states are inde- 
pendently distributed over the whole pure-state space in 
a unitary invariant way? And what is the optimal mean 
probability of success? 

We assume that the input state p is prepared in system 
and N copies of each reference state p a (a = 1,2) in 
systems a\, a>2, ■ ■ ■ , ajv, which will be collectively denoted 
by a. We specify the system which an operator acts on 
by the system number in the parentheses; namely, p(0) 
means that this is an operator acting on system 0, for 
example. 

Our task is then to unambiguously distinguish two 
states Pl (0)pf N (l)pf N (2) and p 2 (0)pf N (l)pf N (2). The 
mean success probability of identification is given by 



/ N Hd) 



2 ^ 

a=l 



tr[E aPa (0)pf N (l)pf N (2)}), (7) 



where {E , Ex, E 2 } is a POVM acting on the whole sys- 
tem and (• • •) represents the average over p\ and pi de- 
fined in the preceding section. When the outcome of the 
POVM is a (= 1, 2), we identify the input p with p a with 
certainty. Outcome of the POVM means we have an 
inconclusive result. Note that the POVM should be in- 
dependent of pi and p2 since we are given no classical 
knowledge on them. 

The average over the reference states can be easily per- 
formed by the use of formula for the average of the n-fold 
tensor product of an identical pure state p 



(p®») = 



(8) 



where S n is the projection operator onto the totally sym- 
metric subs pac e and d n is its dimension, d n — tr [S n ] — 
n+d -iC d -i [13. We find 



pW(d) = 



1 



2djM +1 dN 



tr [EiSN+i (01)5jv(2)] 



-tr [E 2 S N (l)S N+1 (02)} 



(9) 



where <Siv+i(01) is the projector onto the totally sym- 
metric subspace on systems (0,1) = (0, li, la, ■ • • , In) 
and other 5's are defined similarly. 

The POVM should satisfy the following conditions: 

E ,E 1 ,E 2 >0, E + E 1 +E 2 = l (10) 

and, for any p± and p 2 , 

tr [E lP2 (0)pf N (l)pf N (2)] = 0, 

tr[E 2Pl (0)pf N (l)pf N (2)] = 0, (11) 

which implies no error is allowed. It is evident that the 
above no-error conditions are equivalent to 

£iSjv(l)<Siv+i(02) = S N (1)S N+1 (Q2)E 1 = 0, 
E 2 S N (2)S N+1 (01) = 5^(2)5^+1(01)^2 = 0. (12) 
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Now we observe that the set of POVM's satisfying con- 
ditions HI L)| l and l|12l) is convex. Namely, if each of two 
POVM's E a and E' a respects conditions JUJJ and (O, 
so does their convex linear combination qE a + (1 — q)E' a 
for any < q < 1. And the resulting probability, 
Eq.ljnj, is also a convex combination: p(qE+ (1 — g)-®') = 
qp(E) + (1 — q)p(E') in an obvious abbreviated notation. 

We exploit this convexity of the POVM to impose some 
symmetries on the optimal POVM without loss of gener- 
ality. First we notice the problem is symmetric under the 
exchange between systems 1 and 2. Suppose a POVM F a 
is optimal. Then another POVM F' a , defined by 

F[ = TF 2 T, F' 2 = TFiT, Fq = TF T, (13) 

is also legitimate and optimal. Here we introduced the 
exchange operator T between systems 1 and 2. Then a 
new POVM E a = \{F a + F' a ) is also optimal and satis- 
fies the exchange symmetry between systems 1 and 2, 

E 2 = TE\T, E = TE Q T. (14) 

The second symmetry we consider is the unitary sym- 
metry of the distribution of the reference states. If a 
POVM F a is optimal, another POVM defined by 

^(2^+1)^+^(2^+1) (fl = Qj 1; 2) (15) 

is also legitimate and optimal for any unitary operator 
U. Let us construct a POVM by 

E a = J dUU® ( - 2N+1 '>F a (U + )® ( - 2N+1 '> 

(a = 0,1, 2), (16) 

where dll is the normalized positive-invariant measure of 
the group U(d). The new POVM E a is clearly a legiti- 
mate optimal POVM. Furthermore, since E a commutes 
with jj®( 2N + 1 ) for any U, we conclude that E a is a scalar 
with respect to the group U (d) . Thus we can assume that 
the optimal POVM satisfies the exchange symmetry of 
Eqs. (fHE|) and is scalar with respect to the group U{d). 

By the exchange symmetry, the mean probability, 
Eq.@, to be optimized takes the form 

P (N) ( d ) = , 1 , tr[E 1< S JV+1 (01) < S j v(2)] . (17) 

dN+lClN 

And the conditions E\ should satisfy are given by 

Ei > 0, 1 > E x + TE X T, (18) 

and the no-error conditions 

E 1 S N {1)S N+1 {Q2) = S N (l)S N+1 (02)E 1 = 0. (19) 

Finally we note that we can work in the subspace V sym , 
in which each of systems 1 and 2 is both totally sym- 
metric. With this in mind, we set S/v(l) = <5jy(2) = 1 
hereafter. 



IV. CASE OF QUBITS (d = 2) 

In this section we study the case of qubits (d = 2), 
where the individual system can be regarded as a spin- 
1/2 particle and the problem reduces to the angular mo- 
mentum recoupling. In the subspace I4 ym , each system 
a (=1, 2) consisting of N spin- 1/2 particles is totally sym- 
metric, implying the total angular momentum of each 
system is j = N/2. 

We can construct the total angular momentum of the 
whole 2N + 1 systems in two ways. First the com- 
bined system of and 1 has the angular momentum 
Ji = j_ = j - 1/2 or Ji = j + = j + 1/2. Then this 
intermediate angular momentum J\ is coupled with the 
angular momentum j of system 2, resulting in the total 
angular momentum of the whole system J. Using the 
standard notation |14| . we write the resultant eigenstate 
with the total angular momentum J and its z component 
M as 

\A Jl] JM) = \{j^)J l ,r,JM) (J 1 =i_,i+), (20) 

where we ordered three systems as 1 ® (g) 2 on the right- 
hand side. Note that the state | Aj + ; JM ) is totally sym- 
metric in the subspace of systems and 1 and the state 
| Aj_ ; JM ) is not, that is, 

S N+1 (0l)\A j+ ;JM) = \A j+ ;JM), 
S N+ i(01)\Aj_;JM) = 0. (21) 

Another coupling scheme is that systems and 2 are first 
coupled to the intermediate angular momentum J 2 . This 
coupling scheme defines another orthonormal base in the 
whole space, 

\Bj 2 ;JM) = \j,(±j)J 2 ;JM) (J 2 = j+,j-), (22) 

where the three systems are ordered in the same way 
as in Eq. J20J). The state | Bj 2 ; JM ) has the following 
exchange symmetries: 

S N+1 (02)\B j+ ;JM) = \B j+ ;JM), 
Sjv+i(02)| Bj_;JM) = 0. (23) 

For a given set of J(^ 2j + 1/2) and M, the two bases 
(|20|l and l12_?li are related by a unitary matrix, which can 
be taken to be real by the standard phase convention, 

\A Jl ;JM)= Yl Rij 2 \Bj 2 ;JM), (24) 

J2=j+,j- 

where the recoupling coefficient Rj j2 is expressed by the 
Racah coefficient, 



Rj lJ2 = ^(2J 1 + 1)(2J 2 + 1)W(j-Jj; Ji J a ), (25) 
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and its explicit form is given by the following 2x2 or- 
thogonal matrix: 



R J = 



gj+i 

V(2j + J+f)(2j-./+|) 
2j+l 



V(2j + J+|)(2j-J+l) 
2j+l 

2j+l 



(26) 



^{A^JMlEtlAj+iJM) 

JM 



JM 



(31) 



where rows and columns are allocated in the descending 
order of Ji and J2, respectively [l4|. 

Now the no-error conditions 119|) imply that Ei is an 
operator in the space spanned by \ Bj_; JM), which is 
annihilated by 6>jv+i(02) as shown in Eq. (|23|) . Further- 
more, E\ can be assumed to be a L^(2) scalar owing to 
the argument in the preceding section. This means that 
Ei is diagonal with respect to J and is proportional to 
the identity for M. Combining these two properties, we 
find that Ei should have the form 

v-h j 

E i=Y. e -' E \B j _;JM){B j _-JM\, (27) 



where coefficients ej should be non-negative by the pos- 
itivity of E\. 

An upper bound is further imposed on the coefficient 
ej by the remaining condition 1 > E% +TE{T in Eq. lfTHj l. 
This condition can be written as 



2j-| J 



M=-J 



B j _;JM)(B j _;JM\ 
+\A j _;JM)(A j _;JM\) 



v- 



= E e ^ E El Bjl i JM ) °J?J, ( B ^JMl 

J=\ M=-JJ u Ji 

(28) 

where the matrix Ojj is expressed in terms of the re- 
coupling coefficients as follows: 



R 



3-3+ 



R J R J 

\ n 3-3+ n 3-3- 



ID J TDJ 

^3-3+^3-3- 



R 3-3- 



(29) 



Eigenvalues of are readily calculated and found to 
be 1 ± \Rj_j_ |. Thus constraints on the coefficient ej are 
given by 



0<ej < 



1 



1 + \R 



(30) 



3-3- 



Now it is easy to express the trace in Eq. 11171 in terms 
of the recoupling coefficients and ej: 

tr [EiS N+ i{Q\)S N {2)] 



Therefore the probability, Ea. (|17fl . reaches its maximum 
when the coefficients ej takes its upper bound given in 

Thus the optimal mean unambiguous identification 
probability is given by 



2 7+1 2j '~* \ R Li 
v {N) (d = 2) - V V 

^N+1 Z N ~\ 1 + \Rj 



2J+ 1 

2w+i2 



E ''I , )• ( 32 ) 



N 



We used the orthogonality of the recoupling matrix in the 
above derivation. Inserting the explicit form of the re- 
coupling coefficients and performing the sum in the above 
expression, we find a simple formula for pi£2>z(d — 2): 

N 



3{N+ 1) 

The optimal POVM is then given by 



(33) 



23- 



Ei = y - 

E-2 = TEiT, Eq = 1 — Ei — E2 



i— \B j _;JM){B j _;JM\, 

U 3-3- I M=-J 



(34) 



As N goes to infinity, PmaL_(d — 2) approaches 1/3, 
which is equal to the mean unambiguous discrimination 
probability p ma x(d — 2) given in Eq. (J5J. When N = 1, 
on the other hand, the optimal POVM takes the form: 



Ei = \{l - 5 2 (02)), E 2 = H(i - 5 2 (01)), 
Eq = 1 — Ei — E 2 . 



(35) 



which reproduces the one given by Bergou and Hillery 
[Iol | in the case of the equal prior probabilities. For a 
two-spin-l/2-particle system, the state is either symmet- 
ric (triplet state) or antisymmetric (singlet state). There- 
fore, the optimal POVM for N — 1 can also be written 
as 

Ei = ?! $(02) )( $(02) I, E 2 = | $(01)) ($(01) I, (36) 

where we introduced the singlet state | $(01) ) = 
( I ) 1 1 ) — 1 1 ) I ) ) / y/2 for systems and 1 and similarly 
I $(02) ) for systems and 2. 



One might wonder if we really need the complicated 
explicit form i|2ti|) of the Racah coefficients to obtain the 
simple final result of Eq. Ij33|l . Actually we can avoid 
the explicit use of Racah coefficients if we exploit the 
algebraic properties of angular momentum operators. We 
will show it in the general dimensional case treated in the 
next section. 



V. CASE OF ARBITRARY DIMENSION d 

In this section, we generalize the argument in the pre- 
ceding section to the arbitrary dimensional case. The es- 
sential point was the intimate relation between the sym- 
metry properties under system permutations and the an- 
gular momentum of the combined system. The symmetry 
under system permutations is characterized by the rep- 
resentation of the symmetric group S^at+i. And the an- 
gular momentum specifies the representation of SU(2), 
more generally the unitary group U{2). Therefore, in 
the case of arbitrary dimension d, we should classify 
the states according to representations of the symmet- 
ric group S2N+1 and the unitary group U{d). 

Let us introduce the orthonormal base of the total 
space (C d }®( 2N + 1 ) according to irreducible representa- 
tions of the symmetric group S2N+1 and the unitary 
group U(d). We write states in this base as 

I A, a, 6). (37) 

Here A represents an irreducible representation of S^jv+i, 
which is specified by a Young diagram. By the expression 
A = [Ai , A2, • • .] , we denote a Young diagram consisting of 
a set of rows with their lengths given by Ai, A2, . . .. The 
label a indexes orthogonal vectors in a particular S2N+1 
representation space and it runs from 1 to the dimension 
of the 5W+1 representation. It is known that the A also 
specifies irreducible representations of the unitary group 
U(d) and its vectors are indexed by 6, which runs from 1 
to m\{d), the multiplicity of representation A of S2N+1 
on (C d )®( 2N+ V 0. 

As stated before, we can work in the subspace V sym , 
where systems 1 and 2 are both totally symmetric, 
iSat(1) = 1 and Sn(2) = 1- Possible Young diagrams 
A appearing in V sym and the range of the index a associ- 
ated with a particular A can be determined by decompos- 
ing the product of three U (d) irreducible representations 
[1] C§) [N] (g) [N] . We decompose the space V sym into three 
orthogonal subspaces V n (n — 1,2,3) according to the 
number of rows, n, of the Young's diagram (see Fig. 

The subspace V\ consists of totally symmetric states: 

I [2JV + l],b>, b=l,...,m [2N+1] (d), (38) 

where we omitted the index a, since the totally symmetric 
representation of S2N+1 is one dimensional. The states 
in V2 belong to representations specified by Young's dia- 
grams of two rows [Ai, A2], where N + 1 < Ai < 27V and 
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FIG. 1: Decomposition of the product of three U(d) irre- 
ducible representations [1] <X> [N] Cg> [N]. The decomposition 
leads to the three orthogonal subspaces V n (n = 1,2,3) ac- 
cording to the number of rows, n, of the Young's diagram. 

A 2 = 2N + 1 — Ai. Since each of these U(d) represen- 
tations appears twice in V aym , we distinguish the two by 
label a = 1, 2 as follows: 

I [Ai, A 2 ],a,6), a = 1,2, b = 1, . . . , m [Xl ,\ 2 ] (d)< (39) 

The remaining states are those whose Young's diagram 
has three rows and span the subspace V3. 

|[A!,A 2 ,1],6), 6 = l,...,m [All A a ,i](d)» (40) 

where N < X 1 < 2N - 1 and A 2 = 2N - A x . We do not 
need the label a for these states, because each represen- 
tation of this type occurs only once in V sym . Note that 
the length of the third row is always 1. 

Now let us determine a possible form of the POVM 
elements E\ and E%. First of all, E\ should respect the 
no-error conditions, Eq. i|19fl. If \x) is in Vi, it is clear 
that <Sjv+i (02) I x ) — \x). It is also easy to see that 
5jv+i(02)| x) = for I a;) € V3, because representations 
with the Young's diagram of three rows cannot be con- 
structed otherwise. States in V% for a given set of A and b 
can be constructed in two different ways. We can assume 
the label a — 1, 2 for states in V2 is chosen such that 

S N+1 (02)| [Ai,A a ], 1,6) = 0, 

5^+1(02)1 [Ai,A 2 ], 2, 6) = \ [X 1 ,X 2 ],2,b). (41) 

We should also remember that the POVM element E\ 
can be chosen to be a scalar with respect to U(d). All 
these facts lead to the following form for E\ : 

£i=E e l^M] E I A *], hb)([Xi,X 2 ], 1, b I 

Ai b 

+ E e IAi.W]Z)l[^»^,l],6)<[Ai J A2 > l] J 6|» (42) 

Ai b 

where e\ are some non-negative coefficients. The co- 
efficient e\ depends only on A, and the orthonormal 
states I [Ai, A 2 ], 1, b) and | [Ai, A 2 , 1], 6) are complete in 
the space annihilated by the projector Sn+i{02). There- 
fore, this Ei can be most conveniently expressed as the 
following operator form: 

Ei = e(l-«Sjv + i(02)), 

e = E e ^' ( 43 ) 

A 
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where T\ is the projection operator onto the U(d) rep- 
resentation space specified by A. Furthermore, we can 
express E 2 as 



Ej = el-5. 



(01) 



(44) 



by the same operator e, since T\ is symmetric under the 
exchange of systems 1 and 2 and we can assume E% — 
TE\T owing to the conclusion in Sec. II. 

Now that we have determined the possible form of E\ 
and E2, we can proceed to the positivity condition of Eq'. 
namely, 1 > E\ + E%. This condition can be written as 

l>£i+£ 2 = e(2- 5^+1(01) -Sjv+i (02)) 

= e(l-A), (45) 

where we introduced an operator A in the subspace V sym 
to be 



A = 5jv+i(01) +5jv+i(02)-1. 



(46) 



It is convenient to introduce another operator, which 
is the difference of the two projectors: 



D 



5jv+i(01)-5at + i(02). 



(47) 



Note that operators A and D are diagonal with respect 
to A and proportional to identity for the index b, since 
these operators involve just permutation operators. We 
also observe the relations 



A 2 



1 — D 



AD + DA = 0, 



(48) 
(49) 



which can be shown by an explicit calculation using 
5at +1 (01) 2 = 5 w+ i (01) and S N+1 (02) 2 = S N+1 {02). 

The operator A is —1 in the subspace V3, since both 
projectors <Sjv+i(01) and Sn+i(02) annihilate any states 
in V3. In the subspace V2, two eigenvalues of A have 
opposite signs in the invariant subspace associated with 
a given set of A and b. This is because A and D an- 
ticommute and the operation of D changes the sign of 
eigenvalue of A. Note that D does not annihilate any 
state in V2. Combining these facts, we conclude that the 
positivity condition of (|45fl implies the following inequal- 
ity: 



1 + 14 



> e, 



(50) 



in subspaces V% and V3. 

Let us go back to the mean success probability, 
Ea. (|17J) . and write it with the form of E\ given by 
Eq.gSJ. 



v {N \d) = 



1 



IN- 



-id 



N 



-tr [e(l - 5jv+i(02))5jv + i(01))] , (51) 



which can be further rewritten as 
1 



pW(d) = 



tr [e(l-A 2 )}. 



(52) 



In the above equation, we find that the subspaces V\ 
and V2 have no contribution to the trace sum. And |j4| in 
the upper bound of e in Eq. (|50|1 commutes with 1 — A 2 in 
the trace. Therefore, we immediately obtain the optimal 
mean success probability as follows: 



P [N) (d) < 



1 



2dw+idN 
1 

2dN+idN 

(d). 



-tr 



1 + 14 
tr[l-L40 



(53) 



The optimal success probability is thus attained by 
1 



Ex = 



r (l-«5w+i(02)), 



1 + 14 

E-2 = ^-A^(l-<5iV+i(01)), 



E 



1 + 14 
-4 + 14 
1 + 14' 



(54) 



Here we took er^,^,!] = \ l° r simplicity, which is the 
maximum value allowed by Eq. l(5l7|) , though the subspace 

V3 does not contribute to PmaL(d). 

We must still determine eigenvalues of \A\ in the sub- 
space V% in order to evaluate p^sL(d) further. As men- 
tioned above, the operator A is proportional to iden- 
tity with respect to the index b for a given particular 
A. Therefore, eigenvalues of A are independent of the 
dimension d up to multiplicity. Thus we can assume the 
dimension d is equal to 2, which allows us to exploit the 
angular momentum algebra. 

Assuming d = 2, we introduce three sets of angular 
momentum operators: s(0) = |<r(0) for system 0, j(l) 
for system 1, and j(2) for system 2. In the subspace 
V ayia , we have j(a) 2 — ^(t + 1)' since the total angular 
momentum of system a = 1,2 is It is then easy to 
show that the projector Sn+i (0a) can be written in terms 
of the angular momentum operators as follows (a = 1,2): 



Sn+i(0o) 



1 



N + l 



2j(o)-s(0) 



N 



1 



(55) 



Using this form for the projectors, we calculate A 2 . After 
some algebraic calculation involving the Pauli matrices 
and angular momentum commutation relations, we find 



A 2 = 



1 



(N + iy 



(56) 



where J = s(0) + + j(l) is the total angular mo- 
mentum operator. The eigenvalue of J 2 is J(J+ 1) (J = 
I , . . . , N + 1 ) , which implies that eigenvalues of A are 
J+- 

given by ± jfj^ with multiplicity 2 J + 1. 

For a general dimension d, we thus conclude that A in 
V2 has eigenvalues + X N ~^\ with multiplicity Tfi[\ u \ 3 ]((f) , 
since the total angular momentum J in the case of d = 2 
is given by J = = Ai - N - |. 
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P (d) 

max 







□* ++- 
* ++ 



3 .X* 



+x 

X 



xx xx>0<xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx; : 
x d=2 

+ d=3 
* d=4 
□ d=5 



5 10 15 20 25 30 35 40 45 50 



N 



FIG. 2: The optimal mean unambiguous identification proba- 
bility pi^sL(d) as a function of the number of the copies (N) of 

the reference states. As N increases, p^x(d) approaches the 
mean optimal unambiguous discrimination probability shown 
by the horizontal lines. 



Finally we obtain the formula for the optimal success 
probability 



Pmax(^) 



1 



diy + idjy 
2N 



\ 1= N+1 



N+l 



where m^ lt \ 2 ]{d) (A2 = 2N + 1 — Ai) is the multiplicity 
of the S2N+1 irreducible representation [Ai, A2] and given 
by 



(Xi+d- 1)!(A 2 + d - 2)!(Ai - A 2 + 1) 
(d-l)!(d-2)!(Ai + l)!A 2 ! 



(58) 



Let us study the asymptotic value of Pmak(d) when the 
number of the copies N is very large. In this case we can 
replace the sum in Eq.J57|) by a continuous integration 



with respect to x 



N 



1. We find 



(d) -> 2(d-l) / da (1 + xf- 2 {l - x) 
Jo 

2 d ~ 1 (d - IV 
1- ,„ ) (JV-oo), (59) 



(2d- 1) 



ll 



which is equal to p max (d) given by Eq. (J5J. Thus, as 
expected, the unambiguous identification reduces to the 
unambiguous discrimination as the number of the copies 
goes to infinity. Figure [21 displays how the unambiguous 
identification probability approaches the unambiguous 
discrimination probability as the number of the copies 
increases. 

VI. CONCLUDING REMARKS 



We have studied the problem of unambiguously iden- 
tifying the input state of a d-dimensional system with 
one of the two reference states when N copies of each 
reference state are presented with no classical informa- 
tion. We have determined the optimal mean unambigu- 
ous identification probability p[^aL(d) as a function of d 
and N. 

It is interesting to compare the results in this paper 
and those for the identification problem without the no- 
error conditions, which was studied in our previous paper 
0. In both problems the symmetry under system per- 
mutations plays an essential role. This is also true in 
the state comparison studied by Barnett, Chefles, and 
Jex ^l]] j m which one's task is to establish whether or 
not two quantum systems have been prepared in the 
same state. In this case the symmetry under exchanging 
the two systems can characterize the optimal POVM. In 
the state identification, however, we must distinguish the 
two states Pl (0)pf N (l)pf N (2) and p 2 (0)pf N (l)pf N (2). 
Therefore, we must consider the symmetries with re- 
spect to partial permutations among systems and 1 
and among systems and 2. The relevant operators 
are noncommutable projection operators 5/v_|-i(01) and 
<5jv+i(02), which makes the optimization of the success 
probability rather involved. The success probability is 
expressed by the trace of the modulus of some linear com- 
bination of the symmetrizers <S/v+i(01) and <5>jv+i(02): 
D = <5>iv+i(01) — <Sjv+i(02) in the case of the iden- 
tification problem without the no-error conditions and 
A = 5jv+i(01) +5tv + i(02) — 1 in the unambiguous iden- 
tification problem considered in this paper. 

As for the optimal POVM, it was shown that the op- 
timal success probability can be attained by a projec- 
tive measurement in the identification problem without 
the no-error conditions, whereas the optimal POVM ob- 
tained for the unambiguous identification considered here 
is not a projective measurement. 
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